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We present x-ray absorption spectra around the 3d — » 4/ and 4d 4/ excitation thresholds of 
in-plane magnetized Gd and Tb films measured by total electron yield using circularly polarized 
synchrotron radiation. By matching the experimental spectra to tabulated absorption data far 
below and above the thresholds, the imaginary parts of the complex refractive index are determined 
quantitatively. The associated real parts for circularly polarized light propagating nearly parallel 
or antiparallel to the magnetization direction are obtained through the Kramers-Kronig relations. 
The derived magnetooptical parameters are used to calculate soft x-ray reflectivity spectra of a 
magnetized Gd film at the 3d — > 4/ threshold, which are found to compare well with our experimental 
spectra. 
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I. INTRODUCTION 

Magneto-optical (MO) effects in the visible-light re- 
gion are widely used for analyzing magnetic materialsj'^ 
and have important technological applications e.g., in the 
reading process of MO disksi^ They are based on the fact 
that left- and right-hand circularly polarized (CP) light is 
reflected with different intensities from a magnetic mate- 
rial depending on the local magnetization direction. Al- 
though MO effects are small in the visible-light region, 
sensitive detection methods yield enough contrast to dis- 
tinguish bits of opposite magnetization and allow the ob- 
servation of domain structures in optical microscopy:"^ 
One powerful feature of MO techniques is their capa- 
bility to monitor magnetization reversal processes in ap- 
plied magnetic fields, which is not readily achievable with 
electron detection techniques due to Lorentz forces. MO 
techniques in the visible-light region involve optical tran- 
sitions between delocalized valence states, which renders 
it extremely difficult to spectrally separate the magnetic 
contributions of different elements in compounds and 
in advanced layered or nanostructured materials. This 
can be a severe limitation in analyzing magnetic nanos- 
tructures^ or heteromagnetic systems for information 
storage 

Element sensitivity is naturally achieved by employ- 
ing optical transitions that involve core electrons. Large 
magneto-optical signals in the x-ray region were theoret- 
ically predicted for resonant scattering (XRS), i.e., when 
the x-ray photon energy is tuned to the transition energy 
from a core level into a partially filled shell that con- 
tributes to an ordered magnetic moment.^ In fact, XRS 
has been used extensively to study the magnetic struc- 
ture of lanthanide^ and actinide^ materials in the hard 
x-ray regime. 

In the soft x-ray region, large changes in the specu- 
larly refiected x-ray intensity at the ^2,3 edges of tran- 
sition metals upon magnetization reversal have been de- 
tected^Si^^i^ and used in element-specific studies of het- 
eromagnetic systemsi^^d'^d^d^d^d^d^i^" 

In present-day thin films and multilayers, the wave- 



length of soft x rays is comparable to the system dimen- 
sions (typical thicknesses range from 1 to 10 nm). It has 
been shown by several studiesi^ii^iiS that in order to ex- 
tract layer-resolved magnetization profiles from measured 
soft x-ray MO signals of layered systems, a comparison 
with model calculationsSi of reflected intensities is usu- 
ally required. These are based on the Fresnel equations 
and need accurate values of the MO constants as input. 

Several experimental determinations of soft x-ray MO 
constant»ii*22iS^ and reflection coefficients^i have been 
reported for the L2.3 thresholds of ferromagnetic transi- 
tion metals, but results on rare-earth elements have been 
scarce, in spite of the fact that these are often employed 
to achieve high coercive fields in magnetic layers (e.g. 
in spring magnetaS^) or large perpendicular magnetic 
anisotropics.^® Only recently has it been demonstrated 
that sizeable MO signals are obtained from lanthanide 
elements in the soft x-ray region at the N4_5 thresholdsi^i 
In a further study, a huge Faraday rotation has been pre- 
dicted for Gd films at the photon energy corresponding 
to this transition, making use of the experimentally de- 
termined difference in the refractive index for oppositely 
magnetized material)^ 

It is the aim of this paper to perform a quantitative 
determination of the MO constants for the 4d 4/ 
and 3d 4/ absorption thresholds of ferromagnetic 
Gd and Tb metals. Absorption spectra, calibrated by 
matching to tabulated data far from the resonances, pro- 
vide the data basis for a determination of the imagi- 
nary parts of the refractive index. The associated real 
parts are obtained through a Hilbert transformation us- 
ing the Kramers-Kronig (KK) relations. In this sense, 
the magneto-optical parameters derived in this work are 
quantitative; they are consistent with the tabulated val- 
ues far from the resonances, but they have not been de- 
termined in an absolute way. Finally, in order to illus- 
trate the applicability of the MO constants obtained in 
this way, we calculate soft x-ray refiectivity spectra of a 
magnetized Gd film at the M4_5 thresholds and compare 
them with our experimental refiectivity spectra. 
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II. EXPERIMENTAL 

Experiments at the 4c? 4/ absorption thresh- 
olds were performed at the UE56 undulator beam- 
hne^^ of the Berhner Elektronenspeicherring fiir Syn- 
chrotronstrahlung (BESSY II), while those at the 3d — > 
4/ thresholds were performed at the beamlinc ID12- 
B/HELIOS-I of the European Synchrotron Radiation Fa- 
cility (ESRF).^2*2i In the UE56 experiments, the pho- 
ton energy resolution was set to about 100 meV (full 
width at half maximum), which is well below the in- 
trinsic width of the narrow iV4_5 pre-edge absorption 
lines of Gd and Tbi^ By scanning the photon en- 
ergy at slow speed through a synchronized movement 
of monochromator and undulator, an easy normaliza- 
tion of the spectra was made possible; this also allows 
one to exploit the high flux of the undulator beamline 
of about 10i''photons/(s X 100 mA X 0.1% bandwidth) 
over a wide energy range. The degree of circular polariza- 
tion at this Sasaki- type undulator beamline is 98±2%i^ 
The reflectivity measurements at the Gd M4_5 lines were 
also performed at BESSY II. For the absorption measure- 
ments at the M4.5 thresholds at ID12, the energy resolu- 
tion was set to about 0.4 eV and the degree of circular 
polarization was 94±3%.'''^ 
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FIG. 1: (a) Schematic experimental arrangement used for 
measuring absorption of CP light by TEY, using a channel- 
tron, and specular reflectivity by means of a photodiode. In- 
plane magnetized Gd and Tb epitaxial films were prepared on 
W(llO) and measured in situ, (b) XMOKE hysteresis loop 
of an 8 nm-thick Gd film measured at a photon energy of 
147.5 eV. 

Absorption spectra were recorded in total-electron 
yield (TEY) mode using a high-current channeltron. To 
suppress the background of secondary electrons from the 
chamber walls, both the sample and a retarding grid 
placed in front of the channeltron were biased using a 
low- voltage battery. For signal stability, high voltage was 
supplied by a 3.2-kV battery box. The electron-yield cur- 
rent was amplified by an electrometer. TEY spectra of 
the Ad — > 4/ thresholds were normalized to the mirror 
current measured in the last refocusing clement of the 
beamline located in a separate vacuum system in front 
of the experimental chamber. Because it is essentially 
the same physical effect (secondary electron generation) 
that is used for the measurement of both absorption and 
normalization signals, possible distortions due to a lack 
of proportionality between TEY and absorption coeffi- 



cient in the wide photon energy range measured for 
the -/V4^5 transitions are expected to be minimal. The 
film thicknesses were always sufficiently larger than the 
electron's inelastic mean free path (IMFP) relevant for 
TEY measurements so that contributions of the W sub- 
strate to the TEY signal are negligible. In addition, such 
a contribution would be structureless in the measured 
photon-energy ranges and would be removed by the pro- 
cedure of matching the ends of the spectra to tabulated 
values (see below). Refiectivity spectra were recorded 
in specular geometry with a photodiode mounted on a 
rotatable feedthrough; the diode current was measured 
by means of a low-noise electrometer. The experimental 
arrangement is shown schematically in Fig. ^ 

Epitaxial Gd and Tb metal films with thicknesses rang- 
ing from 10 to 50 nm were prepared in situ by vapor de- 
position in ultrahigh vacuum on a W(llO) single-crystal 
substrate. The base pressure in the chamber was in the 
10~^^ mbar range, rising to about 4 x 10"^*^ mbar dur- 
ing film deposition. Compared to transmission methods, 
which give the absorption coefficient in a more direct 
way, TEY detection has the advantage of allowing the 
use of metallic single crystals as substrates. The growth 
of epitaxial films on them is well characterized. In par- 
ticular, controlled annealing of the deposited lanthanide 
films can be performed at the optimum temperatures for 
achieving smooth films with homogeneous thicknesses, 
(for details of film preparation, see Ref. 35). For re- 
manent in-plane sample magnetization, an external field 
was applied along the [110] direction of the substrate us- 
ing a rotatable electromagnet This corresponds to the 
easy axis of magnetization of the Gd and Tb films. In 
the case of the magnetically harder Tb films, the sample 
was cooled from room temperature down to the measure- 
ment temperature (30 K) in an external magnetic field 
on the order of 0.1 T to achieve a single-domain magnetic 
structure with a high remanent magnetization. This has 
been veryfied by means of the magneto-optical Kerr ef- 
fect (MOKE) in the visible-light range in the laboratory 
and in the soft x-ray range (XMOKE) at the beamline. 
As an example. Fig. ^b) shows an element-specific hys- 
teresis loop of a thin Gd film measured by XMOKE. 



III. RESULTS AND DISCUSSION 

A. Absorption 

Figure [21 displays experimental absorption spectra in 
the region of the Gd and Tb iV4_5 (4c; 4/) thresh- 
olds for nearly parallel and antiparallel orientation of the 
magnetization with respect to the spin of the incoming 
CP photons. Figure 13 shows the -M4_5 (3d — > 4/) absorp- 
tion thresholds of both elements. The spectra have been 
corrected for saturation and scaled to fit the tabulated 
values of Henke et alm^ at the ends of the measured pho- 
ton energy ranges following a procedure that is described 
in detail below. 
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FIG. 2: Absorption coefficients of (a) Gd and (b) Tb at tlie 
Ni^5 {Ad 4/) thresliolds of remanentiy magnetized films at 
T = 30 K. CP ligfit was incident ai 9 — 30° witfi respect to tlie 
film plane, i.e. mainly parallel (filled symbols) and antiparal- 
lel (open symbols) to the in-plane sample magnetization. The 
inserts show the regions of the prepeaks measured with higher 
point densities. 
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FIG. 3: Absorption coefficients of (a) Gd and (b) Tb at the 
Af4,5 (3d — + 4/) thresholds of remanentiy magnetized films 
at T = 30 K. CP light was incident at ^ = 30° with respect 
to the film plane, i.e. mainly parallel (filled symbols) and 
antiparallel (open symbols) to the in-plane sample magneti- 
zation. The raw data were taken from Refsi— and—. 




The measured TEY spectra are affected by intrinsic 
saturation, which becomes significant when the atten- 
uation length of the electromagnetic radiation gets 
comparable to de, the inelastic mean free path (IMFP) of 
electrons in the solidi '^^i^^ In that case, the detected sig- 
nal at the channeltron is no longer proportional to the ab- 
sorption coefficient ^. Due to the large absorption cross 
sections at the lanthanide Ad Af and 3d — > 4/ thresh- 
olds, the intrinsic saturation effect is significant and the 
TEY spectra must be corrected for it in order to achieve 
a reliable quantitative estimation of the absorption coef- 
ficient. The saturation correction was made by inverting 
the relations^ 



/io?e -|- sin 9 

for the measured yield Y. Since the constant C depends 
on geometrical parameters of the detection system and 



is not precisely known, we need to determine the value 
of fj,de at some fixed photon energy from independent 
measurements. In the case of Gd M4_5, the saturation 
correction is based on previous results obtained from 
the (6x6) Eu/Gd(0001) surfaceiS It consists of a quasi 
close-packed Eu monolayer on top of the Gd(OOOl) sur- 
face. Since divalent Eu has the same electronic 4/ shell 
configuration (4/^) as trivalent Gd, both elements show 
the same multiplet structure in absorption at the M^^^ 
thresholds. The transitions are well separated in energy. 
A value of (0.30±0.05) for the product (//(ie)pcak of Gd 
at the M5 threshold has been determined by compar- 
ison of the experimental M5 peaks of Gd and Eu; for 
the latter, negligible saturation can be safely assumed 
due to the small thickness of the Eu overlayer (1 atomic 
layer) 4i Furthermore, by comparing the total (energy- 
integrated) M5 absorption signals of Gd and Tb and tak- 
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ing into account the different number of 4/ holes of both 
elements (seven in Gd, six in Tb), we estimate a value of 
(0.20±0.05) for (yLt(ie)peak at the Tb M5 threshold. A con- 
sistency argument can be applied to justify the use of the 
same values for the saturation correction at the 4c? — > 4/ 
thresholds: the relevant electron IMFP at the Ad 4/ 
thresholds corresponds to about 150 eV kinetic energy 
of the Auger electrons which trigger the secondaries de- 
tected at the channeltron; it is typically 3 times smaller 
than at ~ 1200 eV, according to the "universal" curve^^ 
while the values obtained for the absorption length 
change roughly by the same factor (see below). 
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FIG. 4: Optical constants (/3 and 5) for (a) Gd and (b) 
Tb corresponding to the A'4,5 (4ci 4/) thresholds of re- 
manently in-plane magnetized films with the magnetization 
vector nearly parallel (filled symbols) and antiparallel (open 
symbols) to the spin of CP photons. 

The saturation-corrected absorption spectra were 
scaled to fit tabulated values. The photon energy ranges 
of the present spectra are significantly wider than those of 
previous studies (^24^ they include the wide asymmetric 
flanks in the case of the Ad 4/ resonances. This allows 
us to calibrate the absorption spectra by matching both 
ends to the tabulated absorption coefficients^ at photon 
energies where the influence of the resonances is expected 
to be negligible. To this end we fixed the absorption 
coefficients fj,± at the low- and high-energy sides of the 
measured spectra to the values given by the tables of 



Henke et al.^- This procedure defines the ordinate scales 
in Figs. 121 and O The values for the absorption lengths 
/i^^ determined in this way are given in Tables HI and ITU 
The error bars for the quoted values have been carefully 
estimated. They result from an experimental precision 
of ±1% at the considered energy and at both ends of the 
photon energy range, where the spectra were matched 
to the tabulated data. This has been determined from 
the scatter of the data points about their mean values. 
An additional contribution results from the uncertainty 
in the parameter (fide) peak used in the saturation correc- 
tion (see above) . A full error propagation calculation was 
made considering the explicit functional relation of these 
quantities. 



TABLE I: Absorption lengths (in nm) for Gd and Tb at 
the 4d —f 4/ thresholds for nearly parallel and antiparallel ori- 
entations of (in plane) magnetization and photon spin. Values 
are given at the antiresonance (AR), the largest prepeak and 
the maximum of the giant resonance of each element. 
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The absorption spectra at the Ad 4/ thresholds 
display overall shapes of Fano resonances 4^*^ In addi- 
tion, up to 10 eV below the giant resonance, small sharp 
prepeaks can be seen, which correspond to absorption 
final states reached from the respective ground states 
(Gd«57/2,Tb^F6) through violation of the AL = 0, ± 1 
selection rulei^2i2^ 

Note that the x-ray absorption lengths of a few 
nanometers at the Gd and Tb 7V4,5 maxima are the short- 
est in the periodic tablen^i They are accompannied by a 
huge magnetic contrast: the dichroic (difference) signals 
for Gd and Tb amount respectively to 33 and 60 % of 
the nonmagnetic signals, approximated by the average of 
the spectral pairs in Figs. |2| and O 

The absorption cross section determined here for the 
Gd iV4 5 threshold is more than a factor of 2 higher than 
the calculated one for atomic Eu in Ref . \^ These two el- 
ements have the same 4(i^*'4/^ ground state configuration 
and hence the same 4d^4/* multiplet lines, but owing to 
the lower (screened) nuclear charge of Eu one clearly ex- 
pects differences in all parameters of the giant resonance 
spectral profile (height, width and asymmetry). In par- 
ticular the Fano q parameter depends very sensitively on 
the radial matrix elements. Therefore a relative increase 
of the Ad — s- 4/ absorption maximum by a factor of 2 in 
going from Eu to Gd is not unreasonable. In fact a factor 
of 2 is found in other calculations 4^ 

In the case of the lanthanide Af4^5 thresholds {3d Af 
transitions), the matrix elements for Auger decay are re- 
duced compared to the Ad Af resonances. This causes 
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FIG. 5: Optical constants (/3 and 5) for Gd and Tb corresponding to the M4,5 {3d 4/) thresholds of remanently in-plane 
magnetized films with the magnetization vector nearly parallel (filled symbols) and antiparallel (open symbols) to the spin of 
the CP photons. 



TABLE 11: Absorption lengths ^i"^ (in nm) for Gd and Tb 
at the 3d 4/ thresholds for nearly parallel and antiparallel 
orientation of the (in plane) magnetization and the photon 
spin. Values are given at the Ms and M4 maxima of each 
element as well as at the energy equidistant from the two 
maxima: 1200.9 eV for Gd and 1255.6 eV for Tb. 



element magnet. /spin 


Ms 


betw. maxima 


M4 


Gd par. 

antipar. 


8.4 ± 1.0 
14.7 ± 1.7 


280 ± 50 
300 ± 50 


36 ± 4 
27 ± 3 


Tb par. 

antipar. 


6.7 ± 1.3 

15 ± 3 


330 ± 100 
390 ± 130 


33 ± 7 
29 ± 6 



the Fano q parameter to be of the order of 100.'^^ In the 
hmit of large q, the Fano shape approaches a Lorentzian 
and, in fact, the M^^^ absorption line shape is Lorentzian, 
containing hundreds of multiplet components that clus- 
ter into two main groups, the and M4 thresholds'*^ 
(depending on whether the spin of the 3c? hole state is 
oriented paraUel or antiparaUel to the I = 2 orbital an- 



gular momentum). For several characteristic photon en- 
ergies, the values of the absorption length, determined 
after the matching to tabulated data, are presented in 
Table |n] A quantitative determination of these values 
seems particularly pertinent because standard tables^i 
do not include the 3d 4/ transitions, but only the step- 
like 3d — > 6p threshold. Our results for the absorption 
lengths of Gd at the M5 peak are comparable, though 
somewhat larger than the magnetization-averaged value 
measured in transmission in Ref. 49- However, as stated 
by the authors, their absolute values might be affected 
by systematic errors. Note also the huge experimental 
magnetic contrast at M5, which leads to asymmetries of 
40 and 57 % for Gd and Tb, respectively. Some time 
ago, Goedkoop et al. proposed to apply the large mag- 
netization dependence of lanthanide M4_5 absorption for 
constructing line filters to produce CP x rays.^" Further- 
more, the M4.5 lines are clearly separated in energy even 
for neighboring elements in the periodic table like Gd 
and Tb; this represents an advantage for element-specific 
studies. 



6 



B. Magneto-optical constants 

From the values obtained for the absorption coefhcicnt 
we are able to calculate the magneto-optical constants, 
i.e., the real and imaginary parts of the complex index of 
refraction, defined as 



(2) 



where the + and — signs refer to the magnetization point- 
ing either parallel or antiparallel to the CP photon spin 
vector, respectively. The imaginary part is directly re- 
lated to the absorption coefficient through 



1 he 



^^±iE). 



(3) 



The real parts are calculated by means of a Hilbert 
transformation using the KK relations. Due to the bro- 
ken time reversal symmetry inside a magnetized medium, 
the complex refractive index satisfies the symmetry rela- 
tion n±{—E*) = n^(i?), and one has to apply the modi- 
fied KK equations^ 



d+{E)+S^{E) 
S+iE)-S^{E) 



2E 



■ (3+{E') + p^jE') 
^ f3+{E') - P^{E') 



{E' 



E^ 



dE'.{A) 



In order to perform these integrations over the largest 
possible energy range, we resumed to the tabulated val- 
ues for the optical constants outside the measured re- 
gions in the entire range from to 30 keV taken from 
the compilation of Henke et al^ We have included the 
relativistic correction affecting the asymptotic behavior 
of the real part of the atomic scattering factor at large 
photon energies 



f'iE) 



Z* 



Z - (Z/82.5) 



2.37 



(5) 



where Z is the atomic number, as given by the fit of 
Henke et al& to the tabulated values of Kissel and 
Pratti^ The atomic scattering factor and the optical con- 
stants are related through 



5±{E) 
P±{E) 



r^h^e^nat 



2ttE^ 
2ttE' 



f±iE), 

2 mE), 



(6) 



where rg is the classical electron radius and n^t is the 
atomic density of the element. 

The calculated real parts of the refractive index of 
Gd and Tb at the iV4_5 and M4^5 absorption thresh- 
olds are shown in Figs. 0] and [51 for opposite orienta- 
tions of the magnetization and the photon spin vector. 
They exhibit the well-known dispersive behavior, with 
tails ranging far beyond the associated imaginary parts. 
While the latter give the highest magnetic contrast in 



the absorption maxima, the real parts also provide mag- 
netic contrast in regions where the absorption is small. 
This allows the performance of magnetization-dependent 
measurements in reflectivity in regions of different light- 
penetration depths, as the reflected signal is determined 
by both the real and the imaginary parts of the refractive 
index. 

The exact values for n± can only be obtained in an 
ideal experiment where the magnetization direction and 
the photon spin are strictly colinear. Since this is not 
feasible for in-plane magnetized films, there is a contribu- 
tion of transitions with Amj = with a weight l/2sin^ 6 
fRef. ls^ to both the experimentally obtained n+ and n_ . 
This amounts to only a few percent for incidence angles 
up to 30°. Contributions from the opposite magnetiza- 
tion, scaling as 1/2(1 — cos 6')^, can be safely neglected for 
those angles. The contribution of Amj = transitions 
cancels out if only the difference in the MO constants 
— 71- is required, as in, for example, the calculation 
of the Faraday effect 



C. Reflectivity 

As an application of the magneto-optical constants 
determined in this way, we present in Fig. EJa) calcu- 
lated reflectivity spectra of a Gd film in the region of 
the Af4_5 threshold. The reflectivity was calculated using 
the Jones matrix formalism and the Fresnel equations for 
transmision and reflectivity at interfaces of magnetized 
media. ^"^'^^ We considered interference between two chan- 
nels [see insert in Fig.EJa)]: (1) represents the reflection 
at the vacuum/Gd interface and (2) comprises transmis- 
sion through this interface, propagation in the Gd film 
(which includes absorption and Faraday effect), reflection 
at the Gd/W(IIO) surface, propagation back through the 
film, and transmission through the Gd/vacuum inter- 
face. Higher-order paths, including multiple reflections, 
are found to contribute negligibly to the reflected inten- 
sity. In the calculation, we employed the MO constants 
for Gd at the A'h^^ thresholds shown in Fig. [S] together 
with the values for the W substrate taken from Ref. l37i 

Figure |Hfb) shows the corresponding experimental re- 
flectivity spectra of a Gd film on W(llO). The nominal 
film thickness is d = lf±l nm and the light incidence 
angle is 9 = 10 ± 1°. The experimental spectra are well 
reproduced by the calculation shown in Fig.|HJa), includ- 
ing their fine structure. The best agreement between cal- 
culated and experimental spectra was achieved by setting 
d = 12 nm and 6* = 11° in the calculation, in good agree- 
ment with the experimental values. Quantification of the 
experimental reflectivity R was done by normalizing us- 
ing the diode signal in the direct beam. Considering the 
uncertainties in detector position and the simplicity of 
the calculation (for example, no roughness in the film 
was considered) the agreement of the calculated absolute 
intensities with the measured ones can be considered sat- 
isfactory. 
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In conclusion, by quantifying magnetization-depen- 
dent absorption spectra of Gd and Tb, we were able to 
obtain the values of the MO constants of both elements 
at the Ad 4/ and 3d 4/ excitation thresholds in the 
soft x-ray region. For the example of Gd 3d 4/, we 
have shown that the x-ray MO constants obtained here 
can serve as input for reflectivity model calculations of 
thin lanthanide films. In this way, our results open the 
possiblity to apply the power of soft x-ray reflectivity to 
element-specific studies of nanoscaled layered magnetic 
systems containing lanthanide elements. 
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